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SUMMARY 


Funding  to  provide  for  replenishaent  of  reparable  spares  in  the  United 
States  Air  Force  (USAF)  is  appropriated  through  Budget  Progran  1500  (BP-15), 
the  Aircraft  Replenishaent  Spares  Prograa.  BP- 15  includes  the  Peacetiae 
Operating  Stock  (POS)  essential  to  the  peacetiae  readiness  goals  of  the 
Air  Force,  as  well  as  the  War  Reserve  Materiel  (WRM)  needed  to  sustain  forces 
in  a  conflict.  The  requireaents  for  the  POS  and  WRM  portions  of  BP- 15  have 
always  been  coaputed  separately  using  analytical  techniques  that  differ 
widely,  asking  the  interrelationship  between  peacetiae  readiness  and  wartiae 
capability  hard  to  understand  and  quantify  accurately. 

The  Logistics  Manageaent  Institute's  (LMI's)  Aircraft  Availability  Model 
(AAM)  has  been  used  since  1972  by  Headquarters  USAF  in  assessing  the  POS 
requireaaent.  The  AAM  is  a  stochastic,  aulti-echelon,  aulti-indenture  inven¬ 
tory  aodel  that  relates  the  POS  portion  of  BP- 15  to  a  measure  of  materiel 
readiness  called  the  "aircraft  availability  rate."  Consistent  with  its  use  as 
a  long-range  planning  tool  for  peacetime,  it  is  dependent  upon  a  body  of 
"steady-state”  inventory  theory  techniques.  This  working  note  describes  a 
recent  effort  to  extend  the  AAM's  capability  so  that  it  can  assess  aircraft 
availabilities  throughout  a  dynamic  conflict  scenario. 

Because  the  steady-state  mathematics  is  so  deeply  embedded  in  the  exist¬ 
ing  AAM  methodology,  this  working  note  is  necessarily  technical.  It  incorpo¬ 
rates  recent  results  in  the  published  literature  concerning  the  theoretical 
nature  of  resupply  processes  in  a  dynamic  environment,  as  well  as  approxima¬ 
tion  techniques  for  practical  use  of  those  results. 


-  The  aodif  ications  of  the  AAM  described  here  have  been  iaples»nted  as  a 
prototype  in  a  aodel  referred  to  as  the  Surge  aodel.  In  its  current  fora,  the 
Surge  aodel  is  an  assessaent  tool  only;  that  is,  given  a  data  base  of  existing 


stockage  levels,  it  aeasures  the  aircraft  availability  rates  for  a  single 
weapon  systea  throughout  a  specific  scenario  provided  by  the  user.  In  our 
judgaent,  it  represents  a  first  step  toward  a  aore  aabitious  objective  of 
developing  an  integrated  aodel  for  deteraining  the  funding  requireswnts  of 
both  the  POS  and  WRM  portions  of  BP- 15. 
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1 .  BACKGROUND 


The  Logistics  Manageownt  Institute's  (LMI's)  Aircraft  Availability  Model 
(AAM)  relates  the  Peacetiae  Operating  Stock  (POS)  portion  of  Budget 
Prograa  1500  (BP-15),  the  Aircraft  Replenishaent  Spares  Prograa,  to  individual 
weapon  systea  availability  rates. 

Availability  rate  is  defined  as  the  probability  that  a  randoaly  chosen 
aircraft  of  a  particular  aission  design  (MD)  is  not  waiting  for  a  reparable 
spare.  Although  the  AAM  is  noraally  used  in  Prograa  Objectives  Mearoranda 
(POMs)  and  budget  exercises  involving  the  peacetisw  requirement,  this  working 
note  docuaMtnts  a  recent  LMI  effort  to  develop  a  prototype  aodel  for  assessing 
the  perforaance  of  the  supply  systea  under  the  dynaaic  requireaents  of  a 
typical  wartime  scenario.  We  begin  with  an  overview  of  the  structure  of 
BP- 15. 

BP- 15  is  divided  into  two  aajor  categories:  POS  and  War  Reserve  Materiel 
(WRM).  POS  supports  day-to-day  operations,  the  peacetiise  objectives  of  train¬ 
ing,  and  provides  the  materiel  readiness  necessary  in  the  event  of  a  conflict. 
WRM  funding,  on  the  other  hand,  provides  the  additional  stockage  needed  to 
sustain  forces  during  wartiae  levels  of  activity. 

WRM  is  divided  between  prepositioned  and  prestocked  segments.  The  War 
Readiness  Spares  Kits/Base  Level  Self-Sufficiency  Spares  (WRSK/BLSS)  represent 
the  prepositioned  spares  expected  to  cover  the  first  30  days  of  a  conflict. 
WRSK  are  transportable  spares  used  in  support  of  units  to  be  deployed  in  the 
first  30  days;  BLSS  stocks  supplement  the  peacetiae  spares  for  units  that 
engage  in  place.  In  addition  to  the  WRSK/BLSS  prepositioned  segment,  WRM 
includes  the  prestocked  Other  War  Reserve  Materiel  (OWRM)  segment,  which  is 
designed  to  sustain  the  forces  until  the  defense  industrial  base  is  mobilized. 


The  components  of  BP- 15  and  their  roles  are  summarized  in  Figure  1-1. 


FIGURE  1-1.  BP- 15  AND  READINESS/SUSTAINABILITY  OBJECTIVES 
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The  Directorate  of  Logistics  Plans  and  Programs  (LEX)  of  the  Air  Staff  is 
responsible  for  reviewing  the  total  BP- 15  requirement  during  the  Planning, 
Programming,  and  Budgeting  System  (PPBS)  process.  LEX  is  responsible  for  both 
the  preparation  of  the  annual  Budget  Estimate  Submission  (BES)  for  BP-15  and 
the  spares  portion  of  the  Materiel  Readiness  Report  (MRR)  that  is  submitted 
each  year  by  the  Department  of  Defense  (DoD)  to  Congress. 

LEX  manages  the  BP-15  requirement  through  the  Logistics  Capability 
Measurement  System  (LCMS).  The  LCMS  is  an  umbrella  system  that  includes  the 
AAM  as  its  principal  tool  for  analyzing  the  POS  requirement  and  the  Overview 
model  (developed  by  Synergy,  Inc.)  for  analyzing  the  WRM  requirement. 

POS  and  WRM  requirements  are  managed  separately.  Nonetheless,  there  is  a 
clear  relationship  between  existing  peacetime  stocks  and  wartime  capability. 
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Moreover,  there  is  an  operational  sharing  of  POS  and  WRM  stocks  that  affects 
overall  supply  performance.  Spares  are  drawn  from  WRSK/BLSS  stocks  to  support 
the  peacetime  program.  In  fact,  use  of  WRM  spares  in  evaluating  availability 
rates  is  a  choice  open  to  the  AAM  user.  (This  option,  however,  is  normally 
not  exercised  in  setting  POS  requirements . ) 

LMI,  under  the  sponsorship  of  the  Logistics  Concepts  Division  (LEXY)  of 
LEX,  has  developed  a  prototype  model  for  measuring  the  availability  of  a 
specific  weapon  system  that  is  engaged  in  a  given  wartime  scenario.  The 
availability  rate  can  be  calculated  at  any  point  in  time  in  a  given  scenario 
and  is  based  on  existing  inventory  levels  as  recorded  in  the  Recoverable 
Consusqption  Item  Requirements  System  (D041)  data  base.  The  prototype  necessi¬ 
tated  substantial  modifications  of  the  standard,  "steady-state"  methodology 
suitable  in  the  peacetime  setting.  It  also  incorporates  the  effects  of  can¬ 
nibalization,  because  cannibalization  represents  an  acceptable  and  even  essen¬ 
tial  maintenance  policy  in  wartime. 

This  working  note  describes  the  mathematics  of  the  prototype  model  and 
indicates  directions  of  future  development.  It  also  includes  the  major  issues 
involved  in  extending  the  prototype  from  an  assessment  to  a  requirements- 
determination  capability  for  both  POS  and  WRM. 


2.  THE  SURGE  ASSESSMENT  MODEL 


KEY  FEATURES  OF  THE  MODEL 

The  prototype  assessment  model  incorporates  the  key  features  of  the  AAM. 
The  model: 

1.  Considers  the  tradeoffs  between  stocks  at  depot  and  bases.  An 
optimal  base/depot  distribution  is  chosen,  based  on  maximization  of 
the  weapon  system  availability  rate; 

2.  Incorporates  the  levels-of-indenture  relationship  of  the  components. 
Line-replaceable  units  (LRUs)  can  contain  up  to  four  levels  of 
lower-level  shop-replaceable  units  (SRUs).  Shortages  of  LRUs 
directly  affect  availabilities;  SRU  shortages  affect  availabilities 
only  indirectly  by  delaying  the  repair  of  the  higher-level  assem¬ 
blies  ;  and 

3.  Treats  uncertainty  in  the  projected  demands,  condemnations,  etc. 

These  demands  are  portrayed  as  originating  from  a  Poisson  process 
with  an  estimate  for  the  mean  drawn  from  a  Gamma  distribution.  The 
resulting  distribution  for  the  number  in  resupply  has  a  negative 
binomial  form  that  is  completely  specified  by  estimates  of  the  mean 
and  variance-to-mean  ratio  (VMR)  parameters.  A  negative  binomial 

with  a  VMR  value  of  1.0  reduces  to  a  Poisson;  higher  VMR  values 

imply  more  uncertainty  and,  consequently,  less  availability  for  the 
same  level  of  funding.  (See  [9]  for  more  details  of  the  treatment 
of  uncertainty  in  the  AAM.) 

The  Surge  model  is  fundamentally  different  from  the  standard  AAM  in  the 
mathematics  of  calculating  the  various  resupply  pipelines.  (The  term  "pipe¬ 
line"  is  used  throughout  this  working  note  to  denote  the  mean  value  of  the 

number  of  items  in  a  specific  resupply  process.)  The  AAM  uses  the  well-known 
Palm's  Theorem  [10],  which  is  applied  to  every  instance  of  resupply  in  a 
steady-state  situation.  In  addition,  the  multi-echelon  delay  at  a  base 
because  of  depot  backorders  and  the  multi- indenture  delay  of  LRU  repair 
because  of  SRU  backorders  are  dealt  with  through  steady-state  mathematics. 
The  dynamic  nature  of  the  wartime  scenario  calls  for  a  new  approach  that 
incorporates  a  dynamic  version  of  Palm's  Theorem  [4]  and  a  new  perspective  for 


the  expected  backorder  (EBO)  calculation  first  developed  by  Simon  [13].  This 
perspective  results  in  an  algorithm  that  is  theoretically  exact  but  computa¬ 
tionally  difficult.  The  treatment  used  in  the  Surge  model  is  an  extension  of 
the  approximation  technique  developed  by  Slay  [14]  in  a  model  called 
VARI -METRIC. 

In  addition  to  these  fundamental  changes  in  the  nature  of  the  pipeline 
and  EBO  calculations  (which  are  documented  in  Chapter  3) ,  the  Surge  model  has 
been  designed  to  incorporate  the  effects  of  maximum  cannibalization  of  LRUs. 
(It  assumes  no  cannibalization  of  SRUs.)  The  model  evaluates  aircraft  availa¬ 
bilities  for  a  single  weapon  system  with  a  scenario  provided  by  the  user. 

A  typical  scenario  is  shown  in  Figure  2-1.  This  scenario  resembles  those 
in  War  and  Mobilization  Plan  (WMP)  documents  used  by  wartime  planners.  The 
level  of  intensity  usually  appears  in  the  WMP  in  terms  of  sortie-generation 
requirements  per  day.  This  can  be  translated  into  a  scenario  of  flying-hour 
requirements,  given  other  parameter  values,  such  as  sortie  length  and  flying 
hours  per  sortie.  The  approach  taken  here  is  to  assume  a  scenario  that  is 
already  expressed  in  terms  of  time-dependent  flying-hour  requirements;  the 
varying  flying-hour  program  then  leads  to  proportionate  changes  in  failure  and 
condemnation  rates. 

We  have  the  software  developed  for  the  prototype  assessment  model  and 
implemented  it  for  the  C-5A.  The  results  for  this  weapon  system  are  sum¬ 
marized  below. 

C-5A  RESULTS 

The  C-5A  availabilities  are  based  on  the  simple  step-function  scenario 
portrayed  in  Figure  2-2.  The  65  C-5A  aircraft*  are  assumed  to  generate 


The  65  aircraft  represent  the  aircraft  inventory  for  the  second  quarter 
of  1984  as  it  appears  in  PA84-3,  the  Aerospace  Vehicles  and  Flying-Hour 
Programs  System. 


FIGURE  2-1.  TYPICAL  WARTIME  PROFILE 


uniform  demands  at  each  of  three  bases.  The  resulting  plot  of  availability 
against  time  for  the  C-5A  is  shown  in  Figure  2-3.  These  results  are  termed 


"mission  capable"  (MC)  in  that  only  those  items  considered  essential  to  mis¬ 
sion  performance  are  included  in  the  computation. 


FIGURE  2-3.  C-5A  RESULTS 


TIME  (DAYS) 


In  addition,  the  calculation  incorporates  cannibalization  effects  by 
assuming  that  LRUs  are  cannibalized  at  every  base  to  the  maximum  extent  pos¬ 
sible;  no  cannibalization  of  SRUs  during  the  repair  process  is  assumed. 
Figure  2-3  shows  that  C-5A  availability  deteriorates  quickly.  Of  course,  this 
specific  scenario  is  especially  demanding  on  the  relatively  small  number  of 
aircraft  in  the  inventory.  Some  of  the  other  modeling  assumptions  and  pro¬ 
cessing  logic  for  the  C-5A  are  discussed  below  with  a  detailed  treatment  of 
the  mathematical  algorithms  to  follow  in  Chapter  3. 
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The  Surge  prototype  software  operates  on  a  single  mission  design  series 
(MDS).  Some  preliminary  processing  was  required,  therefore,  to  produce  appli¬ 
cation  and  demand  data  for  the  C-5A  alone.  In  particular,  the  C-5A  "share"  of 
existing  assets  for  common  components  (components  applicable  to  multiple 
aircraft  types  including  the  C-5A)  are  prorated  on  the  basis  of  C-5A  usage. 
The  starting  asset  position  (depicted  at  time  T=0  in  Figure  2-2)  is  determined 
in  a  manner  similar  to  the  near-term  availability  calculation  documented  in 
[6] .  All  on-order  assets  are  (optimistically)  included  in  the  asset  position 
calculation  from  which  Figure  2-3  is  derived. 

The  increased  demands  resulting  from  the  surge  are  assumed  to  impact  the 
Organizational  and  Intermediate  Maintenance  (OIM)  program  only.  This  will 
impact  the  depot  with  respect  to  the  number  of  OIM  failures  that  are  not  base 
reparable;  the  Depot  Level  Maintenance  (DLM)  programs,  by  which  we  mean  the 
regularly  scheduled  inspection  and  overhaul  of  operating  components,  is,  in 
the  prototype  model,  assumed  to  continue  at  the  peacetime  rate. 

In  addition  to  weapon  system  availabilities,  the  prototype  software 
identifies  those  specific  components  whose  poor  performance  drives  the 
results.  For  each  component,  it  is  possible  to  compute  the  expected  number  of 
unavailable  aircraft  (ENUN)  resulting  from  inadequacy  in  the  supply  of  spares 
of  that  component. 

Table  2-1  shows  ENUN  values  at  various  times  for  five  of  the  poorest 
performing  components.  Note  that  none  of  these  components  are  problems  in  the 
peacetime  setting.  It  is  also  apparent  that  a  component  ranking  (in  descend¬ 
ing  order  of  ENUN)  depends  on  the  particular  value  of  T.  When  cannibalization 
is  incorporated,  weapon  system  availability  is  driven  by  a  relatively  small 
subset  of  poorly  performing  components.  On  day  30,  for  example,  when  approxi¬ 
mately  62  of  the  65  aircraft  are  expected  to  be  unavailable.  A  single  com¬ 
ponent  accounts  for  57  of  them. 


TABLE  2-1.  PROBLEM  COMPONENTS  FOR  C-5A 
SEPTEMBER  1983  DATA  BASE 


COMPONENT  NSN2 

1 

ENUN  ON  DAY  T 

T=0 

T=20 

— 

T=25 

T=30 

584100 1539566LH 

m 

47.1 

52.7 

57.4 

584100 1453214LH 

ra 

32.2 

40.3 

48.0 

5841000693066LH 

E9 

13.1 

21.4 

32.8 

6610005600303 

7.5 

18.6 

30.7 

6605010182181 

B 

16.6 

22.5 

28.9 

C-5A  TOTAL 

B 

51.8 

57.6 

62.3 

The  ENUN  values  represent  the  expected 
nuaber  of  unavailable  aircraft  expected  to  be 
unavailable  for  lack  of  the  given  coaponent 
(despite  aaxiaua  cannibalization) .  These 
values  are  based  on  a  fleet  of  65  C-5A 
aircraft. 

2 

The  coaponents  listed  are  associated 
with  the  five  highest  ENUN  values  at  day  30. 


The  essence  of  the  Surge  aodel  logic  is  the  calculation  of  the  expected 
nuaber  in  resupply  (pipeline)  at  a  particular  point  in  tiae.  We  use  the 
Dyna-METRIC  approach  described  in  [4]  to  derive  these  values.  The  details  of 
the  pipeline  calculations,  cannibalization  logic,  and  other  aatheaatical 
algorithms  are  described  in  Chapter  3. 
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3.  MATHEMATICAL  ALGORITHMS 


OVERVIEW  OF  THE  COMPUTATIONAL  METHOD 

Because  the  "steady-state"  theory  is  so  deeply  embedded  in  the  AAM's 
computational  methodology,  documenting  the  changes  necessitated  by  the  dynamic 
environment  requires  a  rather  technical  discussion.  We  begin  with  a  brief 
description  of  the  product  formula  for  availability  and  the  assumptions 
implicit  in  that  formula  and  then  discuss  the  incorporation  of  cannibalization 
into  the  availability  calculation. 

With  or  without  cannibalization,  the  availability  calculation  depends  on 
calculation  of  the  backorder  distribution  of  each  component.  This  computation 
is  based  on  the  METRIC  formulation  by  Sherbrooke  [11]  and  Muckstadt's  [8] 
MOD-METRIC  extension  to  handle  levels  of  indenture.  In  METRIC  (and  the  AAM) , 
each  claimant  (base  or  depot)  is  assumed  to  follow  an  (s-1,  s)  inventory 
doctrine.  In  this  context,  the  EBO  calculation  depends  on  the  evaluation  of 
the  inventory  position  s  and  the  probability  distribution  function  (pdf)  for 
the  number  in  resupply. 

In  the  standard  AAM,  we  develop  this  pdf  by  applying  Palm's  Theorem  to 
each  segment  (repair  or  resupply  from  a  higher  echelon)  of  resupply.  Palm's 
Theorem,  a  classical  steady-state  result,  states  that,  given  demands  generated 
by  a  Poisson  process  and  resupply  times  that  are  independent  of  demand,  the 
pdf  of  the  number  in  each  resupply  segswnt  is  itself  Poisson  and  depends  on 
the  mean  resupply  time  but  not  on  the  specific  distribution  of  the  resupply 
process . 

We  include  in  this  chapter  a  detailed  examination  of  the  applicability  of 
Palm's  Theorem  in  a  simplified  but  illustrative  case.  We  then  present  a 


dynamic  version  of  Palm's  Theorem  [1,  4]  and  demonstrate  how  it  is  implemented 
in  the  surge  environment.  Last,  we  summarize  the  computational  modifications 
required  for  the  EBO  calculation  in  a  dynamic  environment. 

THE  PRODUCT  FORMULA  FOR  AVAILABILITY 


We  define  an  aircraft  as  available  if  it  is  not  missing  any  reparable 
coaqponents.  The  availability  rate  for  a  particular  MD  is  the  percentage  of 
the  MD  inventory  that  is  available  at  a  specific  point  in  time.  The  term 
"availability”  is  used  throughout  this  working  note  to  mean  the  probability 
that  an  aircraft  chosen  at  random  from  the  MD  inventory  is  available.  Thus, 
if  an  MD  consists  of  200  aircraft  and  has  an  availability  of  .75: 


Expected  number  of  available  aircraft  =  .75  x  200 

=  150 

Availability  rate  =  150/200  x  100  =  75  percent. 


The  AAM  computes  the  MD  availability  by  means  of  the  product  formula: 


A  =  n.  (1-EB0./(NAC  x  QPA.)), 


where 


EBCh  =  expected  number  of  backorders  for  component  i 
NAC  =  number  of  aircraft  in  MD  inventory. 

QPA^  =  quantity  per  application  of  component  i 

The  product  is  taken  over  all  the  LRUs  installed  on  the  aircraft.  A  more 
thorough  discussion  of  this  formula  can  be  found  in  [9].  Essentially,  we 
derive  it  by  thinking  of  EBO^  as  the  total  number  of  "holes"  (component 
shortages)  that  are  scattered  at  random  across  the  aircraft  within  the  MD 
inventory. 


Equation  (1)  represents  the  evaluation  of  aircraft  availability  without 
consideration  of  cannibalization.  It  assunes  independence  among  the 
individual  component  backorders.  We  shall  now  develop  the  mathematics  for 
incorporating  cannibalization  effects  into  the  availability  calculation. 
CANNIBALIZATION  LOGIC 

The  fact  that  the  AAM  does  not  include  cannibalization  in  its  logic  is 
consistent  with  its  use  as  a  long-range  budget  tool  and  the  widespread  belief 
that  cannibalization,  withdrawal  of  WRM  assets,  and  other  forms  of  expedited 
repair  should  not  be  included  in  the  requirements  determination  process 
itself.  In  measuring  the  ability  of  the  supply  system  to  cope  with  a  specific 
surge  requirement  over  a  relatively  short  period,  however,  it  is  essential 
that  cannibalization  effects  be  included. 

As  described  earlier,  Equation  (1),  for  a  weapon  system  availability,  is 
based  upon  a  random  scattering  of  the  "holes"  (backorders)  across  the  weapon 
system  inventory.  The  effect  of  cannibalization  is  to  consolidate  these  holes 
upon  fewer  aircraft  than  Equation  (1)  implies.  We  adopt  the  approach  of 
assuming  that  cannibalization  is  performed  to  the  maximum  extent  possible  on 
LRUs.  For  a  fixed  point  in  time  and  a  given  MD  under  consideration,  define 
the  random  variables: 

NAVCAN  =  number  of  aircraft  not  available  despite  maximum  cannibali¬ 
zation  of  LRUs 

NAVCAN.  =  number  of  aircraft  not  available  because  of  shortages  of 
component  i . 
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Then  the  cuaulative  distribution  function 
able  NAVCAN  has  the  fora: 


for  the  randoa  vari- 


FNAVCAN(a)  =  Pr<NAVCAN  1  «0  <2> 

=  Pr(MAX.  {NAVCAN.}  <  n) 

=  n.Pr (NAVCAN.  <  n) 

=  n.Pr (BO.  <  n  QPA.), 


where  BO^  and  QPA^  represent  the  nuaiber  of  backorders  and  quantity  per  appli¬ 
cation,  respectively,  for  coaponent  i.  Evaluation  of  Equation  (2)  requires, 
therefore,  soae  knowledge  of  the  probability  distribution  of  the  nuaber  of 
backorders.  Assuaing  this  can  be  done,  the  expected  nuaber  of  aircraft  not 
available  (ENAVCAN)  is  calculated  froa  Equation  (2),  since: 

ENAVCAN  =  ]T  n  Pr (NAVCAN  =  n) 


(1  "  FNAVCAN(n))' 


The  corresponding  calculation  of  availability  with  cannibalization  is  then: 

ACANN  =  (NAC  '  ENAVCAN )/NAC,  (3) 

where 

NAC  =  nuaber  of  aircraft  in  the  weapon  systea  inventory. 

In  suaaary,  the  standard  availability  calculation  without  cannibali¬ 
zation,  as  given  in  Equation  (1),  depends  on  the  expected  nuaber  of  back¬ 
orders  ,  whereas  availability  with  the  cannibalization  foraula,  as  developed  in 
Equations  (2)  and  (3),  depends  on  knowledge  of  the  probability  distribution  of 
backorders.  We  investigate  the  aatheaatics  of  the  backorder  calculation  in 


the  next  section. 


The  presentation  here  does  not  try  to  replicate  the  vast  literature 
concerning  the  METRIC-based  theory  but  does  emphasize  aspects  of  the  method¬ 
ology  that  require  modification  as  we  move  from  steady-state  to  dynamic 

considerations.  We  will  concentrate  especially  on  the  role  of  the  classic 

Palm's  Theorem  [10]  in  derivation  of  the  steady-state  results,  the  dynamic 

version  of  Palm's  Theorem,  as  developed  by  Hillestad  and  Carillo  [4],  and  the 
incorporation  of  VARI-METRIC  [12,  14]  techniques  into  the  AAM  software  in 
order  to  analyze  the  dynamic  effects.  We  begin  with  a  simplified  treatment  of 
the  METRIC  logic. 

We  assiuse  that  the  component  under  consideration  fails  according  to  a 

stationary  Poisson  process.  Notationally ,  we  say  that  the  random  variable  X 
is  Poisson  with  mean  X,  denoted  X  ~  P(X),  if  the  pdf  of  X  has  the  Poisson 
form: 

Pr(X=x)  =  p(x|X)  =  e~k  X*/x!  . 

We  let  p(x|x)  denote  this  Poisson  probability.  The  expected  value  E(X)  and 
variance  VAR(X)  of  X  are  both  equal  to  the  parameter  X.  The  stochastic 
process 

[X^  =  Number  of  occurrences  of  some  event  in  an  interval  of  time  t} 

is  said  to  be  a  stationary  Poisson  process  if: 

1.  The  distribution  of  occurrences  depends  on  the  length  of  the 
interval  and  not  the  end  points; 

2.  Occurrences  in  non-overlapping  time  intervals  of  time  are  independ¬ 
ent  ;  and 

3.  ~  P(At);  i.e.,  Pr(Xt=x)  =  p(xlxt)  =  e  ^t(Xt)X/x!. 

The  parameter  X  is  referred  to  here  as  the  ’’intensity"  or  demand  rate  of  the 
process 
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Suppose  that  we  have  a  component  in  a  multi-echelon  enviromsent  with 
stocks  to  be  distributed  between  a  single  base  and  depot.  Define  the  follow- 
ing  parameters: 

A  =  demand  rate  (per  day) 
r  =  probability  of  repair  at  base 
BRT  =  average  base  repair  time 
OST  =  average  order  and  ship  time 
DRT  =  average  depot  repair  time. 

We  derive  the  pdf  for  the  number  of  backorders  in  the  steady-state  case.  To 
convey  the  essence  of  the  calculation  without  undue  cosq>lexity,  we  assume, 
further,  that: 

1.  All  failures  are  generated  at  the  OIM  or  base  level.  The  base 
failures  may ,  with  a  certain  probability,  be  reparable  at  the  base 
suintenance  activity;  otherwise,  they  are  shipped  to  the  depot  for 
possible  repair.  We  are  not  considering  the  existence  of  a  separate 
DLM  program,  where  components  are  overhauled  on  a  scheduled  basis 
apart  from  the  base  failure  process; 

2.  All  items  shipped  to  the  depot  for  repair  are  in  fact  reparable; 
i.e.,  there  are  no  condesmations;  and 

3.  The  item  has  no  lower-level  subassemblies  (SRUs). 

The  more  general  treatment  is  docuawnted  at  the  end  of  this  chapter. 

In  addition,  we  assume  that  each  claimant  (base  or  depot)  is  stocking  on 
an  (s-1,  s)  operating  doctrine,  where  the  inventory  position  s  defined  by: 

s  =  (quantity  on  hand)  ♦  (quantity  due  in)  -  (quantity  due  out) 

is  kept  constant  by  replenishment,  if  necessary,  on  a  one-for-one  basis.  A 
backorder  exists  if  the  number  in  resupply  (due  in)  exceeds  s.  Specifically, 


II 

» 


the  expected  number  of  backorders  at  any  claimant  operating  under  an  (s-1,  s) 
policy  is  given  by: 

E[B(s)]  =  2  (x-s)p(x),  (4 

x>s 

where 

p(x)  =  the  pdf  for  the  number  in  resupply. 

In  our  simplified  example,  we  develop  the  methodology  for  explicitly 
calculating  the  pdf  p(x)  and  therefore  the  expected  number  of  backorders. 
Moreover,  p(x)  yields  the  pdf  for  backorders  needed  in  the  availability  with 
cannibalization  Equation  (3)  because: 

Pr(l  backorder)  =  p(s+l) 

Pr(2  backorders)  =  p(s+2), 

and,  in  general: 

Pr(k  backorders)  =  p(s+k). 

Consider  first  the  backorders  for  a  random  point  in  time  at  the  depot, 
given  a  depot  inventory  position  s^: 

E^Bd(sd)1  =  X  (x  _  sd)Pd(x)> 


where 


p^(x)  =  Pr(x  units  in  depot  repair). 


Since  base  demands  are  Poisson  with  a  demand  rate  A  and  each  demand  has  a 
probability  l-r  of  being  shipped  to  the  depot  for  repair,  the  demands  at  the 
depot  represent  a  Poisson  process  with  demand  rate  =  A(l-r). 
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We  now  appeal  to  the  classic  version  of  Palm's  Theorem.  Essentially,  it 
says  that  if  demands  are  Poisson  with  intensity  q,  and  the  resupply  process 
has  an  arbitrary  distribution  with  mean  T  with  the  resupply  time  independent 
of  demand,  the  "steady-state"  number  of  items  in  resupply  has  a  distribution 
that  is  itself  Poisson  with  intensity  =  r|T.  (The  independence  assumption 
concerning  resupply  times  and  demand  is,  in  the  queueing  theory  context,  an 
"infinite-server"  assumption.  In  our  setting,  it  is  sometimes  referred  to  as 
the  "slack  repair  capacity"  assumption.)  Thus,  the  pdf  p^Cx)  in  Equation  (5) 
has  the  form: 

p^(x)  =  p(xl\(l-r)  DRT).  (6) 

We  refer  to  the  mean  value  of  a  particular  resupply  process  as  the  pipeline 
value.  Thus,  the  depot  repair  pipeline  DRPIPE  given  by: 

DRPIPE  =  A(l-r)  DRT, 

together  with  the  spares  level  s^  at  the  depot,  completely  specifies  the  pdf 
for  depot  repair,  the  corresponding  pdf  for  depot  backorders,  and  the  expected 
number  of  depot  backorders  given  by  Equation  (5).  We  apply  this  result  in  the 
derivation  of  the  base  backorder  distribution  that  follows. 

Equation  (4)  at  the  base  takes  the  form: 

E[Bb(sb)]  =  2  (x-sb)Pb(x)» 
x>s, 

D 

where 

Pb(x)  =  Pr{x  units  in  resupply  at  the  base}, 

and  sb  is  the  base  inventory  position.  Resupply  to  the  base  has  two  com¬ 
ponents:  (1)  the  quantities  due  in  from  base  repair,  and  (2)  items  due  in 


from  the  depot.  These  resupply  processes  are  independent  because  the  base  and 
depot  demands  are  independent  Poisson  processes  with  intensities  \r  and 
\(l-r),  respectively.  In  any  event,  we  apply  Palm's  Theorem  to  the  entire 
base  resupply  process.  Assuming  that  successive  resupply  times  are 
independent  of  each  other  and  of  demand  (more  on  this  shortly) ,  it  follows 
that: 

Pb(x)  =  p(xlAR), 

where 

R  =  the  mean  base  resupply  time. 

The  problem  reduces  to  finding  R.  The  METRIC  result  for  the  mean  resupply 
time  at  the  base  is: 

R  =  rBRT  +  (1-r) (OST  +  6(sd)DRT),  (8) 

where 

6(sd)  =  E[Bd(sd)]/E[Bd(0)). 

The  reasonableness  of  Equation  (8)  is  easy  to  justify.  With  probability 
r,  a  given  demand  will  be  repaired  at  the  base  in  which  case  its  expected 
resupply  (repair)  time  is  BRT;  otherwise,  it  is  repaired  at  the  depot,  in 
which  case  the  resupply  time  will  be  OST  plus  the  delay  time  at  the  depot 
resulting  from  depot  backorders.  This  depot  delay  clearly  depends  on  the 
stock  level  at  the  depot  and  appears  in  Equation  (8)  in  the  form  6(sd)DRT. 
The  derivation  of  the  6(sd)  equation  above  is  based  on  the  classic  Little's 
formula  of  queueing  theory,  which  states  that,  under  very  general  conditions, 
the  expected  steady-state  queue  length  is  the  product  of  the  demand  rate  and 
the  expected  waiting  time  (see  [3],  for  example). 


Recalling  the  hypothesis  of  Palm's  Theorem,  Equation  (8)  is  also 
dependent  upon  the  assumption  that  the  base  resupply  process  is  independent  of 
demand.  As  stated  earlier,  this  is  the  "slack  repair  capacity"  assumption 
insofar  as  the  repair  processes  at  base  and  depot  are  concerned.  The  depot 
delay  time  is,  however,  dependent  upon  demand  (unless  s^  =  0)  regardless  of 
any  assumptions  about  repair  capacity.  For  example,  the  probability  that  a 
given  demand  is  immediately  filled  from  depot  stocks  (in  which  case  the 
average  base  resupply  time  is  OST)  is  obviously  dependent  upon  the  recent 
demand  pattern  at  the  depot.  The  METRIC  formulation  is  then  an  approximation 
in  that  it  ignores  this  correlation.  The  effect  is  to  understate  backorders. 
(This  issue  becomes  much  more  significant  in  the  dynamic  setting,  where  the 
derivation  just  described  of  the  average  base  resupply  time  breaks  down 
completely. ) 

We  summarize  below  the  procedures  used  in  the  AAM  for  calculating  the 
distribution  of  backorders  for  a  component  under  our  single  base  hypothesis. 
The  AAM,  in  the  standard  implementation,  is  attempting  to  compute  the  weapon 
system  availability  that  corresponds  to  an  expenditure  of  BP- 15  funds  over  a 
specific  period  of  time  referred  to  as  the  "operating  period."  Since  the 
impact  of  monies  spent  is  not  felt  for  a  procurement  leadtime,  the  asset 
positions  and  various  resupply  pipeline  values  for  each  component  are  cal¬ 
culated  for  a  point  in  time  that  represents  an  average  procurement  leadtime 
beyond  the  model  period. 

The  AAM  first  computes  a  total  worldwide  asset  position  (WWASSET)  for 
each  component.  WWASSET  includes  all  assets  projected  to  be  on  hand  or  due  in 
from  various  sources  of  resupply  (base  or  depot  repair,  procurements 
exogeneous  to  the  supply  system,  etc.).  The  AAM  then  finds  the  distribution 
of  assets  between  bases  and  the  depot  that  maximizes  the  availability  of  the 


weapon  system.  Moreover,  it  computes  the  marginal  worth  (as  a  change  in 
availability)  that  results  from  each  additional  procurement  so  that  a  curve  of 
cost  to  availability  can  be  constructed.  In  our  context,  we  are  not  inter¬ 
ested  in  the  marginal  analysis,  only  in  the  evaluation  of  the  backorder  dis¬ 
tribution  for  this  starting  asset  position  WWASSET.  In  our  one  base  example, 

the  AAM  considers  each  combination  (s.  ,  s.)  such  that: 

b  d 


s  +  s  =  WWASSET. 
b  d 


Given 


A  =  (daily)  base  demand  rate 
r  =  probability  of  repair  at  base 
BRT  =  average  base  repair  time 

OST  =  average  order  and  ship  time  (from  depot  to  base) 

DRT  =  average  depot  repair  time, 

we  have  seen  that  the  mean  values  (pipelines)  for  the  various  resupply  pro¬ 
cesses  are  given  by: 

DRPIPE  =  A(l-r)DRT  depot  repair  pipeline 
BRPIPE  =  Ar  BRT  base  repair  pipeline 

OSPIPE  =  A(l-r)OST  order  and  ship  pipeline. 

The  pdf  P^Cx)  for  the  number  in  depot  repair  is,  by  Palm's  Theorem,  Poisson 
with  mean  DRPIPE: 


pd(x)  =  p(xlDRPIPE). 


From  this,  we  can  calculate  the  EBOs  at  the  depot: 


E(Bd(sd)]  =  X  (*-sd)pd(x) 


and  the  depot  delay  function: 

6(sd)  =  E(Bd(sd)]/E[Bd(0)). 

If  we  define  the  depot  delay  pipeline  as: 

DEPDELAY  =  \(l-r)6(s  JDRT, 

a 

we  have,  after  substituting  E{B(sd) ] /E (B(0) ]  for  6 ( sd )  and  noting  that 
E[B(0)]  =  DRPIPE  that  DEPDELAY  =  E[B(sd)].  Palm's  Theorem  applied  to  the 
entire  base  resupply  process  then  yields: 

Pb(x)  =  pb(xlBRPIPE  +  OSPIPE  +  DEPDELAY). 

The  number  of  EBOs  at  the  base  is,  then: 

E[Bb(sb)]  =  X  (x-sb)pb(x). 

x>s, 

b 

Moreover,  the  entire  pdf  for  base  backorders  has  been  determined  so  that  the 
availability  with  or  without  cannibalization  can  be  computed  via  Equations  (3) 
and  (1).  Thus,  the  availability  calculation  depends  upon  the  calculation  of 
the  total  asset  position  WWASSET  and  the  various  pipeline  values:  BRPIPE, 
OSPIPE,  DRPIPE,  and  DEPDELAY.  We  now  discuss  the  analogous  computations  in 
the  dynamic  setting. 

BACKORDER  CALCULATION  IN  THE  SURGE  MODEL  --  SIMPLIFIED  EXAMPLE 


is  a  nonhomogeneous  Poisson  process  with  intensity  function  \(t),  t>0,  if: 

1.  Demands  in  nonoverlapping  intervals  of  time  are  independent;  and 

2.  The  distribution  for  the  number  of  demands  in  the  interval  [tg,tg+t] 

is  Poisson  with  parameter  A(tg+t)  -  A(t^) ,  where  L 

t 

A(t)  =  J*  Mu)  du. 

J  o 

Condition  (2)  can  be  expressed  as: 

Pr{N(tQ+t)  -  N(tQ)  =  x}  =  p(x I A(tg+t)  -  A(tQ))  (9) 

=  exp{-(A(tQ+t)  -  A(tg)) } {A(tQ+t) 

-  A(tQ)]X/x! . 

Thus,  A(t)  represents  the  mean  (and  variance)  for  the  number  of  demands  in  the 
interval  [0,t].  Note  that  if  A(t)  =  \  for  all  t,  A(t)  =  At  and  Expression  (9) 
becomes : 

Pr{N(tQ+t)  -  N(tQ)  =  x}  =  exp{-(Mt0+t)-AtQ)HMt0+t)  *  AtQ]X/x! 

=  exp{-\t}(\t)X/x!  , 

we  have  the  stationary  Poisson  process  described  earlier.  Hillestad  and 
Carillo  [5]  demonstrate  that  Palm's  Theorem  can  be  extended  in  the  following 
sense.  If  demands  occur  according  to  a  nonhomogeneous  process  with  intensity 
function  \(t) ,  and  if  resupply  actions  are  independent  of  each  other  and  of 
the  demand  process,  then  X^.,  the  number  of  items  in  resupply  at  time  t,  is 
itself  a  nonhomogeneous  Poisson  process  with  intensity  p(t)  given  by: 

t 

p(t)  =  J  Mu) [  1-F(u,t) ]  du,  (10) 

o 
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where  F(u,t)  represents  the  cumulative  distribution  function  (cdf)  of  the 
(possibly  nonstationary)  resupply  process;  i.e.: 


F(s,t)  =  Pr{Resupply  action  initiated  at  time  s  will  be  complete 
by  time  t } . 


The  above  expression  for  p(t)  also  assumes  that  the  number  in  resupply  at 
time  t  =  0  is  0,  but  this  is  only  a  notational  convenience.  Note,  however, 
the  key  difference  between  the  dynamic  and  steady-state  versions.  The  average 
number  in  resupply  at  time  t,  p(t),  depends  explicitly ,  in  the  dynamic  case, 
on  the  distribution  of  resupply  times. 

We  return  now  to  the  hypothesis  of  our  simplified  example,  a  component 
with  OIM-generated  failures  only  --  no  SRUs  and  no  condesmations  --  whose 
stock  is  to  be  allocated  between  a  single  base  and  a  depot.  As  before,  let  r 
be  the  probability  of  base  repair  and  assume  that  base  demands  are  nonhomo- 
geneous  Poisson  with  intensity  function  A(t).  Assume  further  that  all 
resupply  times  (BRT,  OST,  and  DRT)  are  deterministic.  The  Dynamic  Palm's 
Theorem  applied  to  the  depot  repair  process  yields  a  depot  repair  pipeline 
DRPIPE(t)  at  time  t,  which,  by  Equation  (10)  is: 


1 

DRPIPE(t)  =  (1-r)  J  [l-FD(u,t))A(u)  du, 


where  F^  is  the  cdf  of  the  depot  repair  process.  But,  under  the  deterministic 
repair  time  assumption: 


FD(u,t)  = 


0  if  t-u  <  DRT 


1  if  t-u  >  DRT, 


v.v.v.- vs 


and  the  preceding  integral  becomes: 


w 

l-r)  J  A(u)  du,  if  t  <  DRT 


DRPIPE(t)  = 


II 

(l-r)  / 


A(u)  du,  if  t  >  DRT. 


t-DRT 


To  siaplify  the  notation,  assume  that  the  aodel  is  in  steady  state  at  tiae 
t  =  0;  i.e.,  deaands  before  the  start  of  the  surge  scenario  are  stationary 
Poisson  with  intensity  A(t)  =  A^  for  t  <  0.  The  steady-state  value  DRPIPE(O) 
is  then: 


DRPIPE(t)  =  Aq( l-r)  DRT,  t  <  0, 


and  Expression  (11)  becoaes ,  for  an  arbitrary  tiae  t: 


DRPIPE(t)  =  /  A(u)  (l-r)  du. 
t-DRT 


The  kth  aoaent  of  the  backorders  at  the  depot  can  then  be  calculated  via: 


E{B*(t)l  =  2  (x-s  Jkp(*lDRPIPE(t)) 


In  particular,  the  expected  value  and  variance  of  depot  backorders  is: 


E(Bd(t)j  =  E(B‘(t)] 


VAR[Bd(t)]  =  E[Bd(t)]  -  (E(Bd(t)])2. 


Unfortunately,  the  steady-state  argtuaent  for  determining  EBOs  at  the  base 
cannot  be  so  easily  extended  to  the  dynamic  case.  There  are  two  reasons: 

-  The  correlation  between  the  depot  demands  and  the  depot  delay  pipe¬ 
lines  is  more  pronounced  than  in  the  steady-state  case;  and 

-  The  previous  derivation  of  the  6(3^)  expression  depended  on  the 
steady-state  Little's  formula.  This0  is  now  wholly  inappropriate. 

The  following  strategy  is  modeled  after  the  approach  taken  by  Simon  [13] 

and  Kruse  [7]  in  developing  the  exact  values  for  EBOs  at  the  base  under 

steady-state  assumptions .  ! 

Let  (X^:t  >  0]  represent  the  stochastic  process  for  the  number  in  base  ! 

resupply  at  time  t.  Then,  because  of  the  deterministic  assumption  for  OST: 

h  *  “t  ♦  DD[t-0ST,t]  +  Bd(t-°ST)'  (16)  I 

where 

i 

f 

I 

BRt  *  number  in  base  repair  at  time  t  j 

DDjt_QST  »  number  of  depot  demands  in  [t-0ST,t] 

B^Ct-OST)  =  number  of  depot  backorders  existing  at  t-OST. 

The  three  random  variables  on  the  right  side  of  Equation  (16)  are  independent  I 

by  the  independent -increment  assumption  in  the  original  demand  process. 

Moreover,  BR^  is  Poisson  with  intensity  A^(t)  given  by: 

\(t)  =  r  Mt), 


and  DD 


[ t-OST, t] 


is  a  Poisson  random  variable  with  parameter: 


Therefore,  the  expected  value  aad  variance  of  BRfc  is  given  by: 


BRPIPE(t)  = 


A.(u)  du. 


while  the  expected  value  and  variance  of  DD ^ t_Qg-p  is: 


DEPDEM(t)  =  (1-r) 


X(u)  du. 


t-OST 


The  pdf  for  depot  backorders  is  known  because: 


Pr{k  depot  backorders  at  time  t-OST} 

=  Prjs^  +  k  items  in  depot  repair  at  time  t-OST} 

=  p(sd  +  k  DRPIPE( t-OST)). 

Thus,  the  pdf  of  is  the  convolution  of  the  three  pdfs  for  BRt>  D1J[t_osT  t]’ 
and  B^Ct-OST).  This  pdf  can  be  computed  exactly  as  in  Kruse  [7],  but  it  is 
difficult  computationally.  Instead,  we  propose  approximating  this  pdf  in  the 
same  way  that  the  VARI-METRIC  formulation  devised  by  Slay  [14}  is  used  to 
improve  the  accuracy  of  the  METRIC-based  EBO  calculation.  Slay  and  Graves  [2] 
have  demonstrated  independently  that  if  the  pdf  for  base  resupply  is  modeled 
as  a  negative  binomial  random  variable,  the  resulting  EBO  calculation  is  much 
improved  over  the  METRIC  approximation.  (The  improvement  can  be  quantified 
because  of  the  exact  expressions  derived  by  Kruse.)  The  pdf  for  a  negative 
binomial  can  be  expressed  in  the  form: 


P(xp  I  ,Q)  = 


where  and  Q  represent  the  mean  and  VMR  of  the  pdf  and  l~(x)  is  the  g 
function  satisfying  the  recurrence  relation: 


r(n+l)  =  nT(n)  =  n! 

for  integer  values  of  n.  Both  Slay  and  Graves  derive  expressions  for  the 
variance  of  the  nuaiber  in  resupply  so  that  the  pdf  can  be  parameterized .  In 
our  context,  the  expected  value  can  be  written,  from  Equation  (16),  as: 

E[Xt]  =  E[BRt]  +  E(DD[t.0ST,t})  +  E [Bd(t-OST ) )  (18) 

=  BRPIPE(t)  ♦  DEPDEM(t)  ♦  E[Bd(t-OST) ] , 

and  the  last  term  is  given  by  Equation  (14).  Similarly,  because  of  the  inde¬ 
pendence  of  the  random  variables  on  the  right-hand  side  of  Equation  (16): 

VAR[Xt]  =  VAR(BRt)  +  VAR(DD[t_osx  t]]  ♦  VAR[Bd(t-OST) ]  (19) 

=  BRPIPE(t)  +  DEPDEM(t)  +  VAR{Bd(t-OST) ) , 

and  the  last  term  is  given  by  Equation  (15).  Then,  p  =  E |  Xt )  and 

Q  =  VAR[Xt]/E(Xt]  are  the  parameter  values  needed  for  the  pdf  Equation  (17). 
We  discuss  briefly  the  extension  of  these  ideas  to  the  sore  general  case 
THE  EBP  CALCULATION  —  EXTENSIONS 

The  preceding  section  details  the  modifications  in  the  pipeline  calcula¬ 
tions  required  in  the  dynamic  environment  so  that  the  EBOs  can  be  computed 
The  arguments  are  presented  for  a  component  stocked  at  one  base,  with  no 
lower-level  assemblies,  no  condemnations,  and  no  DLM-generated  failures  We 
discuss  the  extension  of  our  EBO  methodology  to  the  more  general  case. 

Multiple  Bases 

If  there  are  NB  bases,  NB  >  1,  the  depot  delay  segment  of  the  base 
resupply  process  is  defined  as  the  portion  of  the  depot  backorders  that 


originated  from  the  specific  base.  The  AAM  assuaes  that  all  bases  are  unifora 
with  respect  to  deaand  and  repair  rates;  the  depot  delay  is  therefore  the  saae 
for  all  bases,  and  the  average  base  resupply  tiae  given  by  Equation  (8)  is 
unchanged.  In  the  dynaaic  environment,  the  decoaposition  of  the  base  resupply 
quantity  Xt  given  by  Equation  (16)  is  valid  so  long  as  B^(t)  is  redefined  to 
be : 

B^(t)  =  number  of  depot  backorders  due  to  a  specific  base. 

Simon  (13]  has  shown  that,  for  a  given  number  of  depot  backorders, 
the  backorders  due  to  a  particular  base  are  binomially  distributed  with  proba¬ 
bility  f  equal  to  the  proportion  of  the  total  depot  deaunds  originating  from 
that  base.  The  pdf  for  B^(t)  could  be  computed  explicitly  on  the  basis  of 
this  binomial  splitting.  The  VARI-HETRIC  approach  we  propose  incorporates 
multiple  bases  by  the  following  minor  SKidif ication  of  Equations  (18)  and  (19): 

E ( X  ]  =  BRPIPE(t)  ♦  DEtDEM(t)  ♦  f  E [ B,, ( t-OST)  ]  (20) 

t  <3 

VAR ( X  |  =  BRPIPE(t)  ♦  DEPDEM(t)  ♦  f  (1-f)  E(BJ(t-0ST)]  (21) 

t  a 

♦  f2  VAR ( B ,( t-OST ) ] . 
a 

These  expressions  are  derived  from  the  well-known  formulas  for  the  mean  and 
variance  of  a  binomial  distribution  and  elementary  results  concerning  condi¬ 
tional  expectations  and  variances 
Condeamations 

Condesmations  at  the  depot  can  be  portrayed  as  generating  a 
requirement  from  a  higher  echelon,  i.e.,  the  commercial  vendor.  In  the  steady 
state,  for  example,  suppose,  as  before,  that  a  component  with  a  base  demand 
rate  K  has  probability  1-r  of  being  shipped  to  the  depot  for  repair  Now 


suppose  there  is  a  positive  probability  c  that  the  component  will  be  condemned 


and  ordered1  immediately  froa  the  Manufacturer .  We  have  a  condemnation 
"pipeline"  (scan  value),  which  can  be  written: 

CONPIPE  =  X(l-r)  c  PLT, 


where  PLT  denotes  the  average  procurement  leadtiae.  Equation  (6)  for  depot 
backorders  th^i  takes  the  form: 


ElBd(*d)1  =  X  (x-srf)p(x! DRPIPE  ♦  CONPIPE). 


x>s 


In  the  dynamic  case,  we  get  an  expression  analogous  to  Expression  (13): 

E[B*(t)]  =  X  («-■«,)*  p(*  I  DRPIPE(t)  +  CONPIPE(t) ) , 
x>s. 


(22) 


where 


CONPIPE(t)  =  J  \(u)  ( 1-r)  c  du 
t-PLT 


(23) 


represents  the  expected  condemnations  over  the  interval  [t-PLT, t].  Note  that 
the  base  backorder  calculation  at  time  t  will  depend  on  the  condemnation 
pipeline  value  at  time  t-OST  analogously  to  the  argument  used  in  deriving 
Expression  ( 16) . 

DLM  Demands 

The  Surge  model  now  treats  the  three  depot  level  maintenance  (DLM) 
programs  as  being  unaffected  by  the  dynamic  scenario;  i.e.,  the  DLM  program  is 
assuawd  to  continue  at  the  peacetime  rate. 

^his  (s-1,  s)  assumption  is  not  realistic  at  the  depot.  Actually,  the 
condemnations  should  be  treated  as  reductions  in  the  inventory  position.  This 
distinction  is  insignificant  over  a  scenario  whose  duration  is  substantially 
less  than  the  procurement  leadtime.  The  treatment  here  allows  for  uncertainty 
in  the  condemnations  and  is  also  preferred  for  purposes  of  exposition. 


This  is  appropriate  because  these  programs  are  aot  generated  by  base 
level  activity.  Consider  a  component  with  a  constant  DLM  demand  rate  x.  We 
have  an  additional  depot  resupply  pipeline: 

DLMPIPE  =  x  DRT, 

which  must  be  considered  in  the  depot  resupply  process.  The  backorders  aris¬ 
ing  from  this  DLM  activity  do  not  affect  the  base  backorders  directly.  The 
normal  AAM  processing  logic  handles  the  allocation  of  depot  backorders  to  each 
base.  This  logic  requires  no  modification  in  the  dynamic  environment. 

Levels  of  Indenture 

Muckstadt  [8]  first  developed  a  methodology  for  the  backorder 
calculation  in  a  multi-echelon  and  multi -indenture  inventory  system. 
Muckstadt 's  model,  called  MOD-METRIC,  applied  the  steady-state  Little's 
formula  to  derive  expressions  for  the  additional  increment  in  base  repair  time 
of  a  specific  LRU  due  to  a  backorder  of  a  constituent  subassembly  (SRU). 
Muckstadt' s  approach  was  approximate  in  that  he  ignored  the  correlation 
between  the  base  delay  caused  by  SRU  backorders  and  LRU  demand.  In  any  event, 
the  MOD-METRIC  approach  is  not  applicable  to  the  dynamic  situation,  both 
because  of  this  correlation  effect  and  because  of  the  breakdown  of  the 
Little's  formula  argument.  Instead,  we  apply  a  VARI-METRIC  argument;  i.e.,  we 
develop  expressions  for  the  mean  and  variance  for  the  number  in  base  resupply 
(including,  of  course,  the  effects  of  SRU  backorders)  at  any  time  t  and  model 
the  base  resupply  process  as  a  negative  binomial  distribution  with  these 
parameters.  Unfortunately,  there  is  no  way  of  quantifying  absolutely  the 
improvement  in  this  method,  because  there  is  no  exact  solution  available  in 
the  multi-echelon,  multi-indenture  case.  Sherbrooke  [12]  recently  developed 
the  most  comprehensive  paper  on  the  VARI-METRIC  method  and  demonstrated  its 


success  by  simulation,  but  he  does  not  consider  condemnations  or  non-OIM 
demands  in  his  treatment.  In  the  one-base  case,  the  approach  is  to  decompose 
the  number  in  base  repair  at  time  t,  X^,  as: 

Xfc  =  number  of  LRU  demands  in  {t-BRT,t] 

J 

+  ^  backorders  of  SRU  j  at  time  t. 
j=0 

This  formulation  assumes,  as  does  MOD-METRIC,  that  every  LRU  failure 
corresponds  to  at  most  one  SRU  delay.  If  the  base  repair  time  BRT  is  treated 

as  deterministic,  the  terms  above  are  independent  processes,  since  SRU  back- 

o 

orders  at  time  t  are  the  result  of  LRU  failures  before  time  t-BRT.  In  the 
multi-echelon  treatment,  it  can  be  shown  that  the  number  of  SRU  demands  at  the 
depot  resulting  from  LRU  repairs  at  the  depot  is  binomially  distributed;  from 
this,  it  is  possible  to  derive  expressions  for  the  mean  and  variance  of  the 
overall  resupply  process.  Chapter  4  contains  a  summary  of  the  approach  taken 
in  the  most  general  case. 


2 

Strictly  speaking,  this  argument  assumes  that  the  BRT  for  the  LRU  con¬ 
sists  entirely  of  initial  disassembly  and  fault  isolation  with  negligible  time 
for  reassembly.  In  fact,  this  assumption  also  implies  that  the  SRU  demands 
are  not  "felt"  for  a  length  of  time  equal  to  the  parent  LRU  base  repair  time. 
We  will  return  to  this  point  in  Chapter  4. 


4.  SUMMARY  OF  DYNAMIC  AVAILABILITY  CALCULATION 


We  have  seen  that  the  availability  calculation  depends  fundamentally  on 
the  calculation  of  inventory  positions  (at  each  base  and  the  depot)  and  the 
worldwide  EBO  computation.  The  Surge  prototype  model  estimates  the  total 
asset  position  WWASSET  at  the  asset  cutoff  date  of  the  specific  D041  data  base 
by: 

WWASSET  =  on-hand  serviceable  and  unserviceable  assets 
+  due-in  serviceable  assets 
+  on-order  quantities 
-  due-out  quantities . 

The  presumption  in  this  formulation  is  that  there  is  enough  warning  time  for 

the  unserviceables  to  be  made  ready  by  the  start  of  the  conflict.  The  inclu¬ 

sion  of  on-order  quantities  is  an  option  that  was  exercised  when  the  results 
shown  in  Figure  2-2  were  obtained.  Included  in  the  "on-hand  serviceables"  are 
the  spares  designated  as  WRM  assets. 

Just  as  in  the  AAM,  the  model  considers,  for  a  given  point  in  time,  all 

combinations  of  inventory  levels  at  base  and  depot,  such  that: 

NB 

2  s.  =  WWASSET, 
i=0  1 

where  s  is  the  depot  level  and  s.  is  the  level  at  base  i,  1  <  i  <  NB.  The 

o  i  -  - 

model  chooses  the  combination  that  results  in  the  smallest  number  of  EBOs . 

Involved  in  this  decision  for  a  given  LRU  is  the  allocation  of  stocks  between 
the  bases  and  depot  for  the  constituent  SRUs .  The  tradeoffs  involved  in 

finding  the  stockage  allocations  are  exactly  the  same  as  in  the  normal  AAM 

logic. 


For  given  stockage  levels,  we  summarize  below  the  methodology  for 
computing  the  probability  distributions  for  the  various  resupply  process  in 
the  most  general  case,  i.e.,  in  the  multi-indenture,  multi-echelon  case  with 
multiple  bases,  condemnations,  and  non-OIM  demands  allowed.  The  notation  is 
borrowed  largely  from  Sherbrooke  [12].  We  show  the  results  for  two  levels  of 
indenture  (SRUs  and  LRUs),  but  extensions  to  higher  levels  are  straightfor¬ 
ward.  We  assume  that  the  SRUs  themselves  are  never  delayed  by  shortages  of 
subassemblies,  although  in  practice  this  assumption  need  only  be  made  for  the 
lowest  level  portrayed.  (Five  levels  of  indenture  are  now  considered  in  the 
AAM  and  the  Surge  prototype.) 

Consider  an  LRU  that  is  stocked  at  NB  bases  and  has  J  distinct  SRUs. 
Throughout  this  discussion,  the  index  i  =  0,  1,  2,  .  .  . ,  NB  corresponds  to 
the  possible  "sites'’  for  supply,  with  i  =  0  representing  the  depot  and 
i  -  1,  2,  .  .  . ,  NB  representing  any  of  the  NB  bases.  The  subscript  j 
represents  either  the  LRU  (j  =  0)  or  any  of  the  constituent  SRUs  (j  =  1,  2,  . 

.  .  ,  J).  We  drop  the  uniform  base  assumption;  it  is  not  required  for  this 
development.  Specifically,  we  define  the  following  parameters: 


=  Demand  rate  at  time  t  for  component  j  (j  =  0,  1,  2,  .  .  . ,  J)  at 
site  i  (i  =  0,  1,2,  .  .  .,  NB).  The  base  demand  process  is 
assumed  to  be  nonhomogeneous  Poisson  with  intensity  \„(t). 

=  DLM  demand  for  component  j,  j  =  0,  1,  2,  .  .  .,J. 

=  Probability  that  failure  of  component  j  will  be  repaired  at  site 


Probability  that  LRU  failure  at  base  i  is  caused  by  failure  of 
SRU  j  . 

Deterministic  repair  time  at  site  i  for  component  j,  assuming  no 
delay  for  repair  parts. 

Deterministic  order  and  ship  time  for  component  j  from  depot  to 
any  base. 


to 


p. 

J 


=  Deterministic  procurement  leadtime  for  component  j . 
=  Inventory  position  for  component  j  at  site  i. 


s .  . 

ij 

Some  relationships  hold  among  these  parameters:  The  depot  demand  rate  for  the 
LRU  is  the  sum  of  the  OIM  demands  originating  at  the  bases  and  the  DLM 
demands : 


NB 


X00^  =  ?  \o^t^1"riO^  +  V 

1=1 


while  the  depot  demand  rate  for  the  SRU  demands  resulting  from  LRU  repair  is: 


NB 


KAt) 


Oj 


=  X  A.. 


i=l 


(t) (1-r . . )  +  T,. 


If  X„(t)  =  number  of  component  j  items  in  resupply  for  time  t  at  site  i,  the 


goal  is  to  determine  expressions  for  the  mean  and  variance  of  X„(t),  so  that 


the  parameters  of  the  presumed  negative  binomial  pdf  for  X„  can  be  computed. 


The  EBO  calculation  then  follows. 

We  first  consider  the  resupply  process  for  SRUs  at  the  depot.  Resupply 
in  this  context  is  from  two  sources:  depot  repair  and  replenishment  to  cover 
condemnations.  These  are  independent  Poisson  distributed  processes  at  any 
instant  in  time.  Thus,  for  j  =  1 ,  2 ,  .  .  .  ,  J,  we  apply  the  dynamic  Palm's 
Theorem  to  each  resupply  process  to  get: 


E[Xqj (t) ] 


VAR(XQj(t)]  =  (depot  repair  pipeline) 


+  (condemnation  pipeline) 

Vj(u)r0j  *  du  *  J 


t-R_ 


t-P. 


Vj(u)(1‘roj)du- 


Thus,  the  expected  SRU  backorders  at  the  depot,  E^BCs^)],  can  be  determined. 
Consider  next  the  LRUs  in  depot  resupply  at  time  t,  Xgg(t). 

Xofl(t)  =  (depot  repair  pipeline  for  LRU  at  time  t) 

+  (depot  condemnation  pipeline  for  LRU  at  time  t) 

+  (SRUj  backorders  that  are  delaying  LRU  repair  at  time  t) . 


=  *00(t)  W0 

is  the  fraction  of  SRU  j  depot  demands  that  are  needed  for  LRU  repair.  For 
any  value  B(Sq^)  of  backorders,  Graves  [4]  has  shown  that  the  number  that  are 
delaying  LRU  repair  is  distributed  binomially  with  probability  parameter 
fgj(t).  Since  the  SRU  fault  isolation  is  assumed  to  take  Rgg  days,  the  SRU 
demands  "lag"  the  LRU  demands  by  the  parent  base  repair  time.  We  model  this 
by  looking  at  the  SRU  backorders  at  time  t-RQg.  Combining  this  observation 
with  the  Graves  results,  it  follows  that: 


E(x00(t)1  =  /  \)0(u)r00  du 
t_Roo 


+  /  \)0(u)(1_r00)du 
t-po 


♦  I  tOt(t-R00)Et-R.  JBCs0i>) 


(26) 


VAK[X00(t  >]=/  X00(u)r00du 

t-Roo 


t 


J 

+  ?,  £0j(t-R00)2v",t-R„ 

j=l  J  00  J 


To  obtain  the  mean  and  variance  for  the  SRU  resupply  process  at  the  bases, 
note  that: 

f. .(t)  =  X. . (t) (1-r . .)/\n. 

ij  ij  0j 

represents  the  fraction  of  depot  demand  for  SRU  j  that  is  due  back  to  base  i. 
The  number  of  depot  SRU  backorders  delaying  resupply  to  the  base  is  again 
distributed  binomially  so  that  for  i  =  1,  2,  .  .  . ,  NB  and  j  =  1,  2,  .  .  .,J: 


A..j(u)(l-r.j)du 


(28) 


VAR(X. . (t) ]  =  J  X. . (u)r .  du  +  / 

1J  t-R. .  lJ  lJ  t-T. 


*  VU-VVA*  t.Tl«(.  0j] 


Finally,  we  arrive  at  the  corresponding  expressions  for  the  LRU  resupply 


process  at  the  base.  Let 


fi0  =  Xi0(1“ri0)/X00 


represent  the  fraction  of  depot  LRU  desund  due  back  to  base  i.  Then, 


Elxio(t)1  =  /_  \o(u)riO  du  +  /  XiO(u)(1"riO)dU 

t  “  K .  „  t-T_ 


*  fi0(t-T0)Et-T0IB<*00)' 


VAR[X 


sj 

*  £  *t-«.  „IBUij»- 
j=l  lO  J 

I0(t)1  =  /„  *iO(u)riO  d“  *  /.  \o(u)<1‘r1. 


t-R . 
lO 


+  fi0(t-T0)(1"fi0(t-T0))Et-T0lB(800)l 


+  fi0  (t-T0)VARt-T0lB(B00)l 


+  Z  VAR  R  (B(s..)]. 
j=l  1  *10  lJ 


These  values  are  used  to  paraaeterize  the  pdf  for  X^(t).  The 
expected  LRU  backorders  at  base  i  are  then: 

Et|B(,iO)1  *  I  (x-.i0)Pr(X.0(t)=x), 


and  the  worldwide  backorders  for  the  LRU  at  tiae  t  becoae: 


NB 

EBO  =  Z  E.[B(«.J]. 
i=0 


This  concludes  the  discussion  of  the  EBO  calculation  in  the  dynamic 
environaent. 

LIMITATIONS  AND  ASSUMPTIONS 

Although  the  existing  software  handles  siaple  step-function  scenarios  of 
the  type  used  in  the  C-5A  iapleaentation,  it  can  be  easily  modified  to  handle 
aore  general  scenarios  of  the  type  depicted  in  Figure  2-1. 

The  more  serious  liaitations  of  the  prototype  swdel  are: 

1.  Each  point  plotted  in  Figure  2-2  is  the  result  of  a  separate  Surge 
aodel  run,  with  assets  and  pipeline  values  evaluated  for  that  par¬ 
ticular  day; 

2.  The  aiodel  assumes  that  the  flying-hour  program  requirements  of  the 
scenario  are  always  achieved.  In  the  C-5A  implementation,  there¬ 
fore,  base  failures  are  assumed  to  occur  at  eight  times  the  peace¬ 
time  rate  throughout  the  first  60  days;  in  executing  the  scenario, 
however,  the  nuaber  of  available  aircraft  drops  so  rapidly  that  it 
would  be  iapossible  to  fly  the  desired  program  after  a  few  days;  and 

3.  The  current  version  has  not  incorporated  the  VARI-METRIC  logic;  that 
is,  it  ignores  the  correlation  factors  between  LRU  repair  and  SRU 
backorders  and  between  base  demands  and  depot  backorders.  In  addi¬ 
tion,  it  makes  no  correction  for  the  theoretic  breakdown  of  the 
Little's  formula  arguawnt  for  the  delay  caused  by  backorders  of 
either  higher-echelon  or  lower- indenture  components.  There  is 
reason  to  believe  that  these  errors  are  much  more  significant  than 
in  the  steady-state  case. 


Because  of  Problea  (1),  the  present  software  is  rather  cumbersome  to 
execute.  The  user  must  preselect  the  days  for  which  the  evaluation  is  desired 
and  fill  in  the  rest  of  the  "curve"  by  interpolation  (portrayed  as  linear  in 
Figure  2-2).  This  is  only  an  approximation  and  a  time-consuming  process,  even 
for  a  small  data  base  covering  components  applicable  only  to  the  C-5A. 
Moreover,  there  is  no  linkage  of  the  software  from  one  day  to  the  next.  The 
availabilities  for  each  day  are  calculated  under  the  assumption  of  an  optimal 
distribution  of  spares  between  bases  and  depots  for  that  day.  This  amounts  to 
an  instantaneously  perfect  lateral  resupply  assumption.  The  current  model 
must  be  modified  so  that  the  initial  (day  0)  distribution  of  assets  is  passed 
along  through  time.  It  is  also  not  obvious  how  to  determine  the  initial 
distribution,  because  the  decision  for  peacetime  levels  of  activity  may  be 
different  from  the  wartime  requirement. 

Problem  (2)  can  be  resolved  fairly  quickly,  but  some  parameter  values  are 
required.  Consider  the  C-5A  results.  Figure  2-2  is  based  on  a  peacetime  C-5A 
flying-hour  program  of  4,500  hours  per  month,  which  is  approximately  2.3  hours 
per  plane  per  day.  Under  the  given  scenario,  this  requirement  becomes  18.5 
hours  per  day.  Even  if  aircraft  were  capable  of  flying  24  hours  per  day,  50 
of  the  65  aircraft  would  be  required  in  order  to  fly  the  daily  program;  i.e., 
an  availability  rate  of  at  least  77  percent  is  required.  This  becomes  impos¬ 
sible  by  day  5.  This  problem  can  be  resolved  if  the  maximum  number  of  flying 
hours  per  aircraft  per  day  can  be  supplied.  Of  course,  this  parameter  depends 
on  factors  exogeneous  to  the  supply  system,  such  as  turn  times  or  constraints 
on  crew  or  maintenance  personnel. 

Problem  (3)  must  be  addressed.  Although  VARI-METRIC  techniques  have  been 
coded  into  some  versions  of  the  AAM,  they  have  not  been  thoroughly  tested. 
Incorporating  these  techniques  is  necessary  because  of  the  mathematical  treat¬ 
ment  of  the  EBO  calculation  in  the  dynamic  environment. 


Finally,  the  theory  outlined  in  this  working  note  is  based  on  some 
assumptions.  Perhaps  the  most  limiting  one  is  the  assumption  of  deterministic 
resupply  times.  Sherbrooke  argues  that  only  the  order  and  ship  times  are 
subject  to  this  restriction,  but  it  is  not  clear  that  this  argument  bolds  up 
in  the  dynamic  setting.  This  will  be  an  area  of  future  investigation. 

The  present  SK>del  also  assumes  no  change  1  a  DLM  activity  during  the 
Surge.  We  shall  also  investigate  ways  to  portray  the  wartime  DLM  program  more 
realistically. 


5 .  CONCLUSIONS 


We  have  demonstrated  the  feasibility  of  modifying  the  AAM  so  that  it  can 
evaluate  aircraft  availability  rates  through  time  under  a  given  dynamic 
scenario  of  activity.  The  required  modifications  are  not  trivial,  inasmuch  as 
the  AAM  relies  so  heavily  on  steady-state  inventory  theory  techniques.  The 
algorithms  presented  in  this  working  note  for  assessing  availability  rates  in 
a  dynamic  environment  are  derived  from  a  generalization  of  the  Palm's  Theorem 
result  concerning  the  probability  distribution  of  resupply  processes  and  from 
recent  work  by  Slay  (14),  Graves  (4),  and  Sherbrooke  [12]  in  improving  the 
accuracy  of  the  EBO  calculation  using  VARI-METRIC  techniques. 

We  have  designed  prototype  software  for  executing  this  methodology, 
which,  in  its  current  form,  has  several  deficiencies.  In  continuing  develop¬ 
mental  work  for  the  Air  Staff,  we  plan  to  enhance  and  improve  the  current 
prototype  so  that  it: 

Tracks  the  base/depot  allocation  of  assets  throughout  the  scenario; 

-  Measures  and  adjusts  for  the  infeasibility  of  a  given  scenario  as 
availability  rates  decline;  and 

Incorporates  VARI-METRIC  techniques  crucial  to  accuracy. 

We  also  plan  to  study  ways  in  which  the  underlying  assumptions  on  deter¬ 
ministic  resupply  times  can  be  relaxed  and  DLM  wartime  activity  can  be 
expl  icitly  modeled. 

Finally,  we  emphasize  that,  although  the  work  described  here  is  largely 
evaluative,  we  believe  the  ideas  documented  have  potential  in  a  requirements 
determination  sense,  i.e.,  that  it  is  possible  to  capture  the  functional 
relationship  between  the  WRM  requirement  and  sustainability  in  much  the  same 
way  as  the  AAM  relates  the  POS  requirement  to  readiness  (as  measured  by 


aircraft  availability).  One  of  the  first  questions  related  to  the  WRM 
requireaents  issue  is  an  appropriate  aeasure  of  sustainability.  We  believe 
auch  progress  is  possible  in  this  regard. 

This  effort  represents  a  promising  beginning  to  the  long-range  goal  of 
modeling  in  a  unified  way  the  overall  peacetiae  and  wartime  requireaents  for 
spares . 
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